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^ . Abstract 

^ . 

We present preliminary results of improved measurements of the CP-violating asymmetries and 

D I branching fractions in the decays B^ —>■ tt+tt^, B^ —>■ K^n^, B^ — > vr'^Tr'^, and B^ — > K^ir^ . This 

r^i ' update includes all data taken at the T(4S') resonance by the BABAR experiment at the asymmetric 

PEP-II B-meson factory at SLAC, corresponding to 467 it 5 million BB pairs. We find 
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5^^ = -0.68 ±0.10 ±0.03, 




C^^ = -0.25 ± 0.08 ± 0.02, 




^i^. = -0.107 ±0.016^°;°°^, 




C^o^o = -0.43 ± 0.26 ± 0.05, 


B{B^ - 


-^ 7r°7r°) = (1.83 ± 0.21 ± 0.13) x 10^^ 


3{B^- 


^ i^V°) = (10.1 ± 0.6 ± 0.4) X 10-^ 



where the first error is statistical and the second is systematic. We observe CP violation with a 
significance of 6.7a in B^ —>■ ■k~^'k~ and 6.1(7 in B^ — > K~^-k~ . Constraints on the Unitarity Triangle 



5^ \ angle a are determined from the isospin relation between all B —> -kt^ rates and asymmetries. 
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1 INTRODUCTION 

Large CP-violating effects [1] in the i?-nieson system are among the most remarkable predictions 
of the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing model [2]. These predictions have been 
confirmed in recent years by the BABAR and Belle collaborations, both in the interference of B^ 
decays to CP eigenstates with and without B^-B^ mixing [3-5] and directly, in the interference 
between the decay amplitudes [6] in B^ -^ K^tt^ [5,7]. 

Effective constraints on physics beyond the Standard Model (SM) are provided by high-precision 
measurements of quantities whose SM predictions suffer only small theoretical uncertainties. Both 
experimental and theoretical uncertainties often partially cancel out in the determination of CP- 
violating asymmetries, which makes CP-violation measurements a sensitive probe for effects of 
yet-undiscovered additional interactions and heavy particles that are introduced by extensions to 
the SM. All measurements of CP violation to date are in agreement with the indirect predictions 
from global SM fits [8, 9] that are based on measurements of the magnitudes of the elements of the 
CKM quark-mixing matrix; this strongly constrains [10] the flavor structure of SM extensions. 

The CKM Unitarity Triangle angle a = arg [— ^dKb/KidKibl ^^ measured through interference 
between decays with and without B^-B^ mixing. Multiple measurements of a, with different 
decays, further test the consistency of the CKM model. The time-dependent asymmetry in B^ -^ 
iT~^ir~ is proportional to sin2a in the limit that only the b ^ u ("tree") quark-level amplitude 
contributes to this decay. In the presence of b ^ d ("penguin") amplitudes, the time-dependent 
asymmetry in B^ -^ TT~^7r~ is modified to 

a{At) = [=^ — ^|- p-^ — ^ = S^^ sm (Am^At) - C^^ cos (Am^At) 



|^(At)|2 + |^(At)|2 

^7V7T 



|^|2-|^|2 (1) 



|A|2 + |^|2 



St,j, = v^l - (72^ sin (2q - 2Aa^^) = Vl^-C2^sin2aeflf, 

where At is the difference between the proper decay times of the signal- and tag-side neutral B 
mesons and Arrid is the B^ mixing frequency. Both the phase difference Aotttt = a — acs and the 
direct CP asymmetry Ctttt may differ from zero due to the penguin contribution to the B^ -^ vr+vr" 
decay amplitude A. 

The magnitude and relative phase of the penguin contribution to the asymmetry Stttt may be 
unraveled with an analysis of isospin relations between the B -^ tttt decay amplitudes [11]. The 
amplitudes A^^ of the B -^ vrV-^ decays and A^^ of the B -^ tt^tt^ decays satisfy the relations 



^+0 = -^A+- + AOO, 

\/2 

V2 



(2) 



The shape of the corresponding isospin triangle is determined from measurements of the branching 
fractions and time-integrated CP asymmetries for each of the B -^ tttt decays. No gluonic penguin 
amplitudes are present in the A/ = 3/2 decay B -^ n n^, so, neglecting electroweak (EW) 
penguins, A'^^ = A^^. We define the direct CP asymmetry C^o^^o in B^ -^ tt^tt^ as 

^ \A^r-\A'r (^. 



From the difference in sliape of these triangles for the B and B decay amphtudes, a constraint on 
Aa^TT can be determined with a four-fold ambiguity. 

The phenomenology of the B -^ tttt system has been thoroughly studied in a number of the- 
oretical frameworks and models [12]. Predictions for the relative size and phase of the penguin 
contribution vary considerably, so increasingly precise measurements will help distinguish among 
different theoretical approaches and add to our understanding of hadronic B decays. 

The measured rates and direct CP-violating asymmetries in i? ^ Ktt decays [13-18] reveal 
puzzling features that could indicate significant contributions from EW penguins [19,20]. Various 
methods have been proposed to isolate the Standard Model contribution to this process in order to 
test for signs of new physics. Sum rules derived from {/-spin symmetry relate the rates and asym- 
metries for the decays B^ or B^ to K^tt^, K^tt^, K^tt^, and K^ir^ [21], while SU{3) symmetry 
can be used to make predictions for the Kit system based on hadronic parameters extracted from 
the TTTT system [19]. 

2 THE BABAR DETECTOR AND DATA SET 

The data used in this analysis were collected in 1999-2007 with the BABAR detector at the PEP-II 
asymmetric-energy i?-meson factory at the Stanford Linear Accelerator Center. A total of 467 it 5 
million BB pairs were used. The preliminary results presented here supersede the results in prior 
publications [5, 13, 16]. Roughly 22% more BB pairs have been added to the BABAR data set, and 
improvements have been introduced to the analysis technique, boosting the signal significance. 

In the BABAR detector [22], charged particles are detected and their momenta measured by a 
combination of a five-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) that covers 92% of the solid angle in the T(4S') center-of-mass (cm.) frame, both operating 
in a 1.5-T solenoidal magnetic field. Discrimination between charged pions, kaons, and protons 
is provided by a combination of an internally reflecting ring-imaging Cherenkov detector (DIRC), 
which covers 84% of the cm. solid angle in the central region of the BABAR detector and has a 91% 
reconstruction efficiency for pions and kaons with momenta above 1.5 GeV/c, and the ionization 
(dE/dx) measurements in the DCH. Neutral-cluster (photon) positions and energies are measured 
with an electromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) crystals. The photon energy 
resolution is cte/E = |2.3/-E(GeV)^'^ © 1.9} %, and the angular resolution from the interaction 



point is ae = 3.9°/^/E{GeY). 

3 ANALYSIS METHOD 

Many elements of the measurements discussed in this paper are common to the decay modes 
B^ -^ h+h''{h = TT ov K), B^ ^ tt^tt^, and B^ -^ K°7r°. The signal S-meson candidates (B^ec) 
are formed by combining two particles, either tracks or vr" or Kg candidates. The event selection 
differs for each mode, and is described in detail below. 

The number of B decays and the corresponding CP asymmetries are determined in extended 
unbinned maximum likelihood (M.L.) fits to variables described below. The likelihood is given by 
the expression 

C = exp 




■ M 



(4) 



where the product is over the number of events N, the sums are over the event categories M, rii is 
the coefficient for each category as described below, and the probabihty-density function (PDF) V 
describes the distribution of the variables x in terms of parameters a. The PDF functional forms 
are discussed in Sec. 3.3.1, 3.3.2, and 3.4. 

3.1 Track and K° Selection 

For particle identification in the B^ -^ h^h'^ sample, we make use of the track's Cherenkov 
radiation in the DIRC as well as its ionization energy loss dE/dx in the DCH. 

For the DIRC information to be used, we require that each track have the associated Cherenkov 
angle (^c) measured with at least six signal photons detected in the DIRC, where the value of Oq is 
required to be within 4.0 standard deviations from either the pion or kaon hypothesis, which effec- 
tively removes any candidate containing high-momentum protons. Electrons are explicitly removed 
based primarily on a comparison of the track momentum and the associated energy deposition in 
the EMC, with additional information provided by DCH dE/dx and DIRC 9c measurements. 

The ionization energy loss dE/dx in the DCH is used either in combination with DIRC informa- 
tion or alone, which enables a 35% increase in the B^ -^ h'^h'~ reconstruction efficiency compared 
to using only the tracks with good DIRC information. A detailed DCH dE/dx calibration that 
we developed for the B^ -^ h~^h'~ analysis takes into account variations in the mean value and 
resolution of dE/dx values with respect to changes in the DCH running conditions over time and 
the track's charge, polar and azimuthal angles, and number of ionization samples. The calibration 
is performed with large high-purity samples (> 10® events) of protons from A -^ Ptt^, pions and 
kaons from D*+ -^ D^tt^ {D^ -^ K^tt^), and Kg -^ vr+vr" decays that occur in the vicinity of the 
interaction region. 

Kg -^ TT+vr^ candidates are reconstructed from pairs of oppositely charged tracks. The two- 
track combinations are required to form a vertex with a x^ probability greater than 0.001 and a 
7r+7r~ invariant mass within 11.2 MeV/c^ (3.7cr) of the Kg mass [23]. 

3.2 7r° Selection 

We form vr" — > 77 candidates from pairs of clusters in the EMC that are isolated from any charged 
tracks. Clusters are required to have a transverse energy deposition consistent with that of a photon 
and to have an energy E^ > 30 MeV for B° -^ 7r°7r° and E^ > 50 MeV for B° -^ -K's^°- We use 7r° 
candidates that fall within the invariant-mass range 110 < rriyy < 160 MeV/cP. 

For the B^ -^ tt^tt^ sample, we also use tt^ candidates from a single EMC cluster containing 
two adjacent photons (a merged n^), or one EMC cluster and two tracks from a photon conversion 
to an e^e~ pair inside the detector. To reduce the background from random photon combinations, 
the angle 9^ between the photon momentum vector in the vr'^ rest frame and the t:^ momentum 
vector in the laboratory frame is required to satisfy | cos^^j < 0.95. The tt^ candidates are fitted 
kinematically with their mass constrained to the nominal t:^ mass [23]. 

Photon conversions are selected from pairs of oppositely charged tracks with an invariant mass 
below 30 MeV/cP whose combined momentum vector points straight away from the beam spot. The 
conversion point is required to lie inside the detector material. Converted photons are combined 
with photons from single EMC clusters to form tt^ candidates. 

Single EMC clusters containing two photons are selected with the transverse second moment, 
S = ^^ Ei X (Aui)'^ /E, where Ei is the energy in each CsI(Tl) crystal and Aa^ is the angle 



10 



between the cluster centroid and the crystal. The second moment is used to distinguish merged vr" 
candidates from both single photons and neutral hadrons. 

3.3 Event Selection in B^ -^ tt+tt-, B° -^ K+tz', and B^ -^ Tr^Tr^ 

Two kinematic variables are used in the B^ -^ h^h'" and B^ -^ tt^'k^ analyses to separate B- 
meson decays from the large e^e~ -^ qq {q = u, d, s, c) combinatoric background [22]. One is 

the beam-energy-substituted mass m-ES = \/('5/2 + Pj • pbY /Ef — p^, where ^/s is the total e^e~ 
cm. energy, (£'j,pj) is the four-momentum of the initial e^e~ system in the laboratory frame, and 
ps is the laboratory momentum of the B candidate. The other is /S.E = E'^ — \/s/2, where E^ is 
the B candidate's energy in the cm. frame. 

Two additional quantities take advantage of the event topology to further separate B decays 
from the qq background. The absolute value of the cosine of the angle 9s between the sphericity 
axes [24] of the B candidate's decay products and that of the remaining tracks and neutral clusters 
in the event, computed in the cm. frame, is peaked at 1.0 for the jet-like qq events but has a 
flat distribution for B decays. We require |cos^s| < 0.7 for B^ -^ tt^tt^ and [cos^sl < 0.91 
for B^ -^ h'^h'~. For the B^ -^ h'^h'~ sample, we further require that the second Fox- Wolfram 
moment [26] satisfy R2 < 0.7 to remove a small remaining background from e~^e~ -^ t~^t~ events. 

To improve the discrimination against qq events, a Fisher discriminant T is formed as a linear 
combination of the sums Lq = ^^ |p*| and L2 = ^j |p*| cos'^6*, where p* are the momenta and 
6* are the angles with respect to the thrust axis [25] of the B candidate, both in the cm. frame, 
of all tracks and clusters not used to reconstruct the signal S-meson candidate. In the case of 
B^ -^ TT^TT^, we improve the sensitivity of the signal by combining J^ with three other quantities in 
a neural network. These are the | cos0g| described above, | cos9b\, where 9b is the angle between 
the center-of-mass momentum vector of the signal B and the beam axis, and | cos^rl, where 9t is 
the angle between the thrust axis of the signal S-meson's daughters and the beam axis. 

3.3.1 S° -^ TT+TT- and B° -^ K+tt' 

We reconstruct candidate decays Srcc -^ h^h'^ from pairs of oppositely charged tracks in the 
polar-angle range 0.35 < ^lab < 2.40 that are consistent with originating from a common decay 
point with a x^ probability of at least 0.001. The remaining particles are examined to infer whether 
the other B meson in the event (-Btag) decayed as a B^ or B^ (flavor tag). We perform an unbinned 
extended M.L. fit to separate B^ — > tt+tt" and B^ -^ K'^it~ decays and determine simultaneously 
their CP- violating asymmetries S'tttt, Ctttt, and Ak-k and the signal and background yields and 
PDF parameters. The fit uses particle-identification, kinematic, event-shape, i?tag fiavor, and At 
information. 

The variables t^es and /S.E are calculated assuming that both tracks are charged pions. The 
B^ -^ TT+TT^ events are described by a Gaussian distribution for both w-es and /S.E, where the 
resolutions are found to be 2.6 MeV/c^ and 29 MeV, respectively. For each kaon in the final state, 
the IS.E peak position is shifted from zero by an amount that depends on the kaon momentum, with 
an average shift of —45 MeV. We require mES > 5.20 GeV/c^ and |Ai?| < 0.150 GeV. The large 
region below the signal in rriEs effectively determines the background shape parameters, while the 
wide range in /^E allows us to separate B^ decays to all four final states (tt+tt", X+tt", tt+X", 
and K~^K~) in a single fit. 

We construct 9c PDFs for the pion and kaon hypotheses, and dE/dx PDFs for the pion, kaon. 
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Figure 1: The average difference between tlie expected values of DIRC 9c and DCH dE/dx for pions 
and kaons at 0.35 < 0iab < 2.40, divided by the uncertainty, as a function of laboratory momentum in 
B° -^ K+TT^ decays in BaBAR. 



Table 1: Average tagging efficiency e, average mistag fraction w, mistag fraction difference Aw — w{B'^) — 
w{B'^), and effective tagging efficiency Q for signal events in each tagging category. The quantities are 
measured in the large-statistics -Bflav sample of fully reconstructed neutral _B-meson decays. 



Category 


e{%) 


w{%) 


Aw (%) 


Q(%) 


Lepton 


8.96 ± 0.07 


2.9 ± 0.3 


0.2 ± 0.5 


7.95 ± 0.11 


Kaon I 


10.81 ± 0.07 


5.3 ± 0.3 


0.0 ± 0.6 


8.64 ± 0.14 


Kaon II 


17.18 ± 0.09 


14.5 ± 0.3 


0.4 ± 0.6 


8.64 ± 0.17 


KaonPion 


13.67 ± 0.08 


23.3 ± 0.4 


-0.6 ± 0.7 


3.91 ± 0.12 


Pion 


14.19 ± 0.08 


32.6 ± 0.4 


5.1 ± 0.7 


1.73 ± 0.09 


Other 


9.55 ± 0.07 


41.5 ± 0.5 


3.8 ± 0.8 


0.28 ± 0.04 


Total 








31.1 ±0.3 



and proton hypotheses, separately for each charge. The K~tt separations provided by Oq and dE/dx 
are complementary: for Oq, the separation varies from 2.5a at 4.5 GeV/c to 13a at 1.5 GeV/c, while 
for dE/dx it varies from less than l.Ocr at 1.5 GeV/c to 1.9a at 4.5 GeV/c (Fig. 1). For more details, 
please see Ref. [5]. 

We use a multivariate technique [27] to determine the flavor of -Btag- Separate neural networks 
are trained to identify leptons and kaons from B and D decays and soft pions from D* decays. 
Events are assigned to one of seven mutually exclusive tagging categories (including untagged 
events) based on the estimated average mistag probability and the source of the tagging information. 
The quality of tagging is expressed in terms of the effective efficiency Q = ^^ efc(l — 2wk)'^., where 
Ek and Wk are the efficiencies and mistag probabilities, respectively, for events tagged in category 
k. The difference in mistag probabilities is given by Aw = w^o — w^o- Table 1 summarizes the 
tagging performance measured in a large data sample of fully reconstructed neutral -Bfiav decays to 
DW-(7r+,p+, a+). 
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The time difference At = /S.z/ (3^c is obtained from tlie known boost of the e^e~ system 
(/37 = 0.56) and the measured distance Az along the beam (z) axis between the Srec and i?tag 
decay vertices. A description of the inclusive reconstruction of the i?tag vertex is given in [28]. We 
require |At| < 20 ps and o"At < 2.5 ps, where UAt is the error on At determined separately for each 
event. The signal At PDF for B^ -^ 'k^it~ is given by 



4±(At„,eas) = " ,_" <! {IT^W) 



e-|At|/r 
"4r 

± (1 - 2wk) [S^^ sin (Am^At) - C^„ cos (Am^At)] } ® R{At^c^ - At), 

(5) 

where /^ ifk) indicates a B^ (B^) flavor tag and the index k indicates the tagging category. 
The resolution function i?(Atmcas — At) for signal candidates is a sum of three Gaussian functions, 
identical to the one described in Ref. [28], with parameters determined from a fit to the -Bflav sample 
(including events in all seven tagging categories). The background At distribution is also modeled 
as the sum of three Gaussians, where the common parameters used to describe the background 
shape for all tagging categories are determined simultaneously with the CP parameters in the 
maximum likelihood fit. 

The M.L. fit includes 28 components: B^ signal decays and background with the final states 
TT^vr^, K^TT^ , K^TT^, and K^K^ where either the positively charged or the negatively charged 
track, or both, have good DIRC information (2 x 4 x 3 = 24 components) plus the p7r~, pK^ , vr+p 
and K~^p background components where the (anti)proton has no DIRC information. The K^tt^ 
event yields are parameterized as n^±^^ = uk-k (1 T -^^kT) /^- ^^^ other coefficients are products 
of the fraction of events in each tagging category, taken from i^fiav events, and the event yield. 
The background PDFs are a threshold function [29] for ttt-es and a second-order polynomial for 
AE. The J- PDF is a sum of two asymmetric Gaussians for both the signal and background. We 
used large samples of simulated B decays to investigate the effects of backgrounds from other B 
decays on the determination of the CP-violating asymmetries in B^ -^ it~^it~ and B^ -^ K^tt~ 
and determined them to be negligible. 

3.3.2 B° -^ ttOttO 

B^ -^ TT^TT^ events are identified with an M.L. fit to the variables w-es, AE, and NN, the output of 
the event-shape neural network. We require niES > 5.20 GeV/c^ and \AE\ < 0.2 GeV. Tails in the 
EMC response produce a correlation between rriES and AE, so a two-dimensional PDF, derived 
from detailed Monte Carlo (MC) simulation, is used to describe signal. The NN distribution is 
binned in ten bins (equally populated for signal) and described by a parametric step-function PDF 
with 9 height parameters taken from the MC and fixed in the fit. -Bflav data are used to verify 
that the MC accurately reproduces the NN distribution. The qq background PDFs are a threshold 
function [29] for rriESj a second-order polynomial for AE, and a parametric step function for NN. 
For qq events, A^A'^ is not distributed uniformly across the bins but rises sharply toward the highest 
bins. We see a small linear correlation between the shape parameter of the to-es threshold function 
and the AA bin number, and this linear relation is taken into account in the fit. All qq background 
PDF parameters are allowed to float in the M.L. fit. 

The decays B+ -^ p+7r° and 5° -^ A:°7r° {K^ -^ 7r°7r°) add 71 ± 10 background events to 
B^ -^ TT^vr^ and are included as an additional fixed component in the M.L. fit. We model these B- 
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decay backgrounds with a two-dimensional PDF to describe rriEs ^^'^ AE, and with a step function 
for NN, ah taken from MC simulation. 

The time-integrated CP asymmetry is measured by the S-flavor tagging algorithm described 
previously. The fraction of events in each tagging category is also constrained to the corresponding 
fraction determined from MC simulation. The PDF coefficient for the B^ -^ tt^tt^ signal is given 
by the expression 

^fkKo^o[l - s,(l - 2x)(l - 2wk)C^o^o}, (6) 

where fk is the fraction of events in the tagging category k, N^o^o is the number of B^ -^ tt^tt^ 
decays, x = 0.188 ±0.003 [23] is the time-integrated B^ mixing probability, and Sj = -|-1(— 1) when 
the Stag is a B^ (^O). 

3.4 B° -^ K°7rO 

For each B^ -^ i^^vr'^candidate, two independent kinematic variables are computed. The first 
one is mB, the invariant mass of the reconstructed B meson, i?rcc- The second one is mmissj 
the invariant mass of the other B, Stag, computed from the known beam energy, by applying 
a mass constraint to i?rcc [30]. For signal decays, tub and m^iss peak near the B^ mass with 
resolutions of ~ 36 MeV/c^ and ~ 5.3 MeV/c^, respectively. Both the w-miss and tub distributions 
exhibit a low-side tail due to the leakage of energy deposits out of the EMC. We select candidates 
within the ranges 5.11 < remiss < 5.31 GeV/c^ and 5.13 < tub < 5.43 GeV/c^, which include a 
signal peak and a "sideband" region for background characterization. In the events with more 
than one reconstructed candidate (0.8% of the total), we select the candidate with the smallest 
X^ = X]i=7r0 ii"" ("^j ~ "^D^/^mi) where rrii {m[) is the measured (nominal) mass and ami i^ ^^^ 
estimated uncertainty on the measured mass of particle i. 

We exploit topological observables, computed in the cm. frame, to discriminate jet-like e'^e" -^ 
qq events {q = u, d, s, c) from the nearly spherical BB events. In order to reduce the number of 
background events, we require L2/L0 < 0.55, where Lj = YlilPil'^^^'^ ^i ^^^ 9* are computed 
with respect to the sphericity axis [24] of the Srec candidate. We compute cos0^, the cosine of 
the angle between the direction of the B meson and the nominal direction of the magnetic field {z 
axis). This variable is distributed as 1 — cos^ 9'^ for signal events and is nearly flat for background 
events. We select events with | cos^^j < 0.9. We also use the distributions of L2/LQ and cos 6*^ to 
discriminate the signal from the residual background in a M.L. fit. Using a full detector simulation, 
we estimate that our selection retains (34.2 it 1.2)% of the signal events, where the error includes 
both statistical and systematic contributions. The selected sample of B^ -^ K^tt^ candidates is 
dominated by random K^tt^ combinations from e^e~ -^ qq fragmentation. Using large samples of 
simulated BB events, we find that backgrounds from other S-meson decays are small, 0(0.1%); 
however, we study in detail the effect of a number of specific B decay channels. The dominant ones 
are B^ -^ p'^K^, B^ -^ K*^tt^, and B^ -^ Kgir^T:^, and we include this effect in our study of the 
systematic errors. 

For the B^ -^ K^n^ decay, where no charged particles are present at the decay vertex, we 
compute the decay point of the -Brec using the knowledge of the Kg trajectory from the measurement 
of vr+vr^ momenta and the knowledge of the average interaction point [31]. 

We extract the signal yield from an extended unbinned M.L. fit to kinematic, event-shape, 
flavor-tag, and decay-time quantities. The use of tagging and decay-time information in the M.L. 
fit further improves discrimination between signal and background. We have verified that all 
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correlations are negligible, and so construct the likelihood function as a product of one-dimensional 
PDFs. Residual correlations are taken into account in the systematic uncertainty, as explained 
below. 

The PDFs for signal events are parameterized based on a large sample of fully reconstructed 
B decays in data and from simulated events. For background PDFs, we select the functional form 
from the background-dominated sideband regions in the data. 

The likelihood function is defined as: 

C(Sf,Cf, Nsig , iVbkg , fsig , /bkg , a) = ^^-j (7) 



X 



n [^sig/sig4g'^sig(^»> y«; Sf, Cf) + iVbkg/bkgeCkg'^bkg(^j, y*; «)] 



ieg 



n [iVsig(l - fsisK,^VU^{x^■ Cf) + 7Vbkg(l - /bkg)e&kgnkg(2^«; a)] 



X 
ie6 



where the A^ selected events are partitioned into two subsets: i £ g events have At information, 
while i G b events do not. Here, /sig (/bkg) is the fraction of signal (background) events G g, 
and 1 — /sig (/bkg) is the fraction of events G b. The probabilities Psig and Pbkg are products 
of PDFs for signal (sig) and background (bkg) hypotheses evaluated for the measurements Xi = 
{tub, mmiss, L2/L0, cos 9%, flavor tag, tagging category} and jji = {At, a At}- "Pgig and P^j^g are 
the corresponding probabilities for events without At information. In the formula, a represents the 
set of parameters that define the shape of the PDFs. Along with the CP asymmetries Sf and Cf, 
the fit extracts the yields Ngig and A^bkg) the fraction of events fsig and /bkg, and the parameters a 
that describe the background PDFs. 

4 RESULTS AND SYSTEMATIC UNCERTAINTIES 

4.1 B^ ^ 7r°7r° Results 

Results from the M.L. fit for the B^ -^ tt^tt^ decay mode are summarized in Table 2. Distributions 
of rriESj A.E, and NN for B^ -^ vr'^Tr^ are shown in Fig. 2, where a weighting and background- 
subtraction technique, s'Plots [32], is used to display the signal events. The same technique is used 
to display the qq background as well, shown in the insets. 

The uncertainty in the efficiency for the B^ -^ tt^tt" decay mode is dominated by a 3% sys- 
tematic uncertainty per vr'', estimated from a study of r ^ tttt^Ut decays. There is an additional 
1.0% uncertainty in the resolution of the signal shape and a 0.45% uncertainty due to the limited 
knowlegde of the w-es and AE peak positions in data, estimated by shifting the rriES and AE 



Table 2: Results for the S'^ -^ tt'^tt'^ and B''^ -^ /C^tt'^ decay modes: signal yields A^sig, efficiencies, branching 
fractions, and time-integrated CP asymmetries. When two uncertainties are given, the first is statistical and 
the second systematic. 

A'^gig Efficiency Branching fraction Asymmetry 

B^ -^ vr'^Tr^ 247 ± 29 (28.8 ±1.8)% (1.83 ± 0.21 ± 0.13) x 10"*^ -0.43 ± 0.26 ± 0.05 
B^ ^ K^TT^ 556 ±32 (34.2 ± 1.2)% (10.1 ± 0.6 ± 0.4) x 10"^ [33] 
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Figure 2: gVlots for B^ -^ tt^tt^ signal (background shown in the inset plots): (top left) m-Es, (top right) 
AE, (bottom) the binned AW. The line in each plot shows the corresponding PDF. 



means and resolutions by amounts determined from MC-data comparison in a control sample of 
B~^ -^ 7r"'"7r'' events. We also take an uncertainty of 1.5%, determined from the -Bflav sample, due 
to the I cos^sl requirement. Systematic uncertainties involving the M.L. fit are evaluated by vary- 
ing the PDF parameters and refitting the data. These contribute an uncertainty of 8.3 events to 
the branching-fraction measurement and an uncertainty of 0.05 to C^Ot^o. The various systematics 
sources are tabulated in Table 3. 



4.2 B° 



TT+TT" 



and B° -^ K+tt- Results 



Results for the B^ -^ h^h'^ decay modes are listed in Table 4. The correlation coefficient between 
5,r7r and Ctttt is found to be —0.056, and the correlation between C-j^n and Ar-k is 0.019. In Fig. 3, 
we show s'Plots for ttt-es, Ai?, and J- for the B^ — > h'^h'~ signal and background. The direct CP 
asymmetry in B^ -^ K~^tt~ is apparent in the distribution of AE plotted separately for B^ and B^ 
decays, shown in Fig. 4. We show the distributions of At for B^ — > K n^ signal and background 
decays in Fig. 5. In Fig. 6, we show the distribution of At separately for B^ -^ ■k'^tt~ events tagged 
as B^ or i?^, and the asymmetry a(At). The central values and errors for S^^j^ and C-n-,^ are shown in 
Fig. 7, along with confidence-level contours corresponding to statistical significances ranging from 
1 to 7 standard deviations. Our measurement excludes the absence of CP violation in B^ -^ tt^tt^ 

into account). 



0) at a confidence level of 2 x 10 , or 6.7a (where systematic errors are taken 
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Table 3: Systematic uncertainties in the determination of the B° -^ tt^tt^ signal yield (A''^o^o) and branch- 
ing fraction, and the direct CP asymmetry C^o^o . The total branching-fraction systematic is the sum in 
quadrature of the uncertainties on the signal yield, the signal efficiency, and the i3-meson counting. 



Source 


iV^o^o 


<ysyst{B)/B 


CttOttO 


Peaking background 
Tagging 

Background shape 
Signal shape 


±4.9 
±0.35 

±5.5 
±3.8 




±0.030 
±0.034 
±0.023 
±0.020 


Total fit systcmatics 


±8.3 


3.4% 


±0.055 


7r° efficiency 
cos 9 s selection 
neutrals resolution 
Toes and AE shape 


6.0% 
1.5% 
1.0% 
0.5% 


Total efficiency systcmatics 


6.3% 


Number of BB pairs 


1.1% 


Total systematic error 




7.2% 


±0.055 



Table 4: Results for the B" -^ h^h' decay modes. For each mode, the number of signal events A^sig and 
CP asymmetries are shown. Statistical, followed by systematic, uncertainties arc given for the asymmetries. 



Mode 



B^ -^ vr+vr 



K 



Slg 



Asymmetry 



TT 



1394 ±54 S"^^ = -0.68 ±0.10 ±0.03; C^^ 
5410 ±91 ^i^^ = -0.107 ±0.016 



0.25 ±0.08 ±0.02 

0.006 
-0.004 



Table 5: Summary of absolute systematic errors on Aktt- The total is calculated as the quadrature sum of 
each contribution. To address the Aktt bias due to hadronic interactions of charged kaons with the detector 
material, we shift the Akti value obtained in the fit by +0.0050. 



Source 



Uncertainty 



Material interactions 

Be and dE/dx PDFs 

Potential MC bias 

Alternative DIRC parameterization 



±0.0053 -0.0025 
0.0020 
0.0011 
0.0016 



Total 



-0.0060 -0.0037 



17 




5.26 5.28 5.3 
m„s (GeV/c^) 



0.1 
AE (GeV) 



2 

Fisher discriminant 



Figure 3: The distributions of (left) niEs, (middle) AE, and (right) Fisher discriminant T: (top) 
background-subtracted for S" — > tt+tt^ signal, (middle) background-subtracted for i?° -^ K^tt~ signal, 
(bottom,) signal-subtracted for all h^h'^ background candidates in the data. The curves represent the PDFs 
used in the fit and reflect the fit result. The structure to the left of the signal AE peak for i?° -^ tt+tt^ is 
consistent with the expected background from other charmless modes, which is negligible above —0.10 GeV. 
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Figure 4: The background-subtracted distribution of AE for signal K^n^ events, comparing (solid) B^ 
and (dashed) B^ decays. 
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Figure 5: (Left) the background-subtracted distribution of At for signal K^tt^ and (right) the signal- 
subtracted At distribution for background candidates in the data. The curves represent the PDFs used in 
the fit and reflect the fit result. 



Table 6: Summary of systematic uncertainties on S^^r and C^r 



Source 


Jtttt 


^TTTT 


DIRC 9c 


0.0064 


0.0050 


DCH dE/dx 


0.0032 


0.0037 


Signal At 


0.0199 


0.0055 


SVT local alignment 


0.0004 


0.0002 


Boost/z scale 


0.0021 


0.0013 


PEP-II beam spot 


0.0028 


0.0014 


B flavor tagging 


0.0146 


0.0138 


Amd, tbo [23] 


0.0004 


0.0017 


Potential bias 


0.0041 


0.0043 


Doubly Cabibbo-suppressed decays [35] 


0.007 


O.OfG 


Total 


0.027 


0.023 
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Figure 6: The background-subtracted distributions of At for signal vr^vr" events tagged as (top) 
B^ or (middle) B^, and (bottom) their asymmetry a{At) (Eq. 1). The curves represent the PDFs 
used in the fit and reflect the fit result. 
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Figure 7: S'tttt and Ctttt in B^ -^ vr^vr : the central values, errors, and confidence-level (C.L.) 
contours for 1 - C.L. = 0.317 (Ict), 4.55 x lO^^ (2^7), 2.70 x lO'^ (3ct), 6.33 x IQ-^ (4cr), 5.73 x lO^^ 
(5(t), 1.97 X 10~^ (60") and 2.56 x 10~^^ (^o"), calculated from the square root of the change in the 
value of — 21n£ compared with its value at the minimum. The systematic errors are included. The 
measured value is 6.7(7 from the point of no CP violation (S-,^-,^ = 0, C^tt = 0). 
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Systematic uncertainties for the direct CP asymmetry Ak-k are listed in Table 5. Here, Ak^w is 
the fitted value of the K^tt^ event-yield asymmetry A^^ shifted by +0.0051q'qq| to account for a 
bias that arises from the difference between the cross sections of K^ and K~ hadronic interactions 
within the BABAR detector. We determine this bias from a detailed MC simulation based on 
GEANT4 [34] version 7.1; it is independently verified with a calculation based on the known 
material composition of the BABAR detector [22] and the cross sections and material properties 
tabulated in Ref. [23]. The corrected K^tt^ event-yield asymmetry in the background, where 
no observable CP violation is expected, is consistent with zero: —0.005 it 0.004 (stat)l^Q'gQ3 (syst). 
Systematic uncertainties for the CP asymmetries S'tttt and Ctttt are listed in Table 6. They are 
dominated by uncertainties in the parameterization of B-flavor tagging and vertexing, and (for 
Ctttt) in the effect of CP violation in Stag [35]. 

4.3 B° -^ K°7r° Results 

Results for the B^ -^ K^tt^ decay mode are summarized in Table 2. In Fig. 8, we show gVlots 
for rriiniss, m^, L2/L0, and | cos^ljl for signal events, with background distributions shown in the 
insets. 

To compute the systematic error associated with the statistical precision on the parameters of 
the likelihood function, we shift each parameter by its associated uncertainty and repeat the fit. 
For At and the tagging parameters, the uncertainty is obtained from the fit to the -Bflav sample, 
while for the other parameters it is obtained from MC; the total error is obtained by summing the 
individual contributions in quadrature. This fit systematic also accounts for the limited statistics 
available to determine the shape of the likelihood function in Eq. 7. We find a systematic error 
of 1.2 events on the K^tt^ yield. As an additional systematic error associated with the data-MC 
agreement of the shape of the signal PDFs, we also quote the largest deviation observed when the 
parameters of the individual signal PDFs for remiss, m^, L2/L0, and cosO*^ are floated in the fit. 
This gives a systematic error on the yield of 2.5 events. The output values of the PDF parameters 
are also used to assign a systematic error to the selection efficiency of the cuts on the likelihood 
variables. Comparing the efficiency in data to that in the MC, we obtain a relative systematic 
error of 1.5%. We do not assign a systematic uncertainty on the scale of rrijaiss and ms because we 
float these variables in the fit. We evaluate the systematic error due to the neglected correlations 
among fit variables using a set of MC experiments, in which we embed signal events from a full 
detector simulation with events generated from the background PDFs. Since the shifts are small 
and only marginally significant, we use the average shift in the yield (-1-2.2 events) as the associated 
systematic uncertainty. 

We estimate the background from other B decays to be small in the nominal fit. We account for 
a systematic shift induced on the signal yield by this neglected component by embedding simulated 
B background events in the data set and evaluating the average shift in the fit result: -1-5.2 events 
on the signal yield. We adjust the signal yield accordingly and use half of the shift as a systematic 
uncertainty. 

For the branching fraction, additional systematic errors come from the uncertainty in the selec- 
tion efficiency, the counting of BB pairs in the data sample (1.1%), and the branching fractions in 
the B^ decay chain, B{Kl -^ vr+vr-) = 0.6920 ± 0.0005 and B{t:^ -^ 77) = 0.98798 ± 0.00032 [23]. 
The systematic uncertainties are summarized in Table 7. 
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Figure 8: Distributions of (a) mmiss, (b) niB, (c) LiJLq, (d) cos^J for background-subtracted events in 
the B^ —^ K'^TT^ sample. The solid curves represent the shapes of signal PDFs as obtained from the M.L. 
fit. The insets show the distributions and PDFs for signal-subtracted data. 



Table 7: Summary of dominant contributions to the systematic error on the measurement oiB{B'^ — > K'^tt' 
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5 CONCLUSIONS 

The CP-asymmetry and branching-fraction results described in this paper are: 

S^^ = -0.68 ±0.10 ±0.03, 
C^^ = -0.25 ± 0.08 ± 0.02, 
^x. = -0.107 ±0.016+°:°°^, 
C^o^o = -0.43 ± 0.26 ± 0.05, 
B{B^ -^ 7r°7r°) = (1.83 ± 0.21 ± 0.13) x 10"^ 
B{B^ -^ K^-K^) = (10.1 ± 0.6 ± 0.4) X 10"^ 

We combine B{B^ -^ ti^tt^) with the branching fractions B{B^ -^ vr+vr-) = (5.5 ±0.4 ±0.3) x 10"^ 
and B{B^ -^ 7r=*=7rO) = (5.02 ±0.46 ±0.29) x 10"^ previously measured by BABm [13,14] to evaluate 
the constraints on both the penguin contribution to a and on the CKM angle a itself. Constraints 
are evaluated by scanning the parameters of interest, |Aa7r7r| = |a— ctcffl and a, and then calculating 
the x^ for the five amplitudes {A^^ ^ A^ , AP^ , A^ , A^^) from our measurements and the isospin- 
triangle relations [8]. The x^ is converted to a confidence level (C.L.) as shown in Fig. 9. The 
upper bound on | ActTTTrl is 43° at the 90% C.L., and the range [23°, 67°] in a is excluded at the 90% 
C.L. If we consider only the solution preferred in the SM [37], a is in the range [71°, 109°] at the 
68% C.L. Somewhat more restrictive new constraints on a have been found in the measurements 
ofB^pp and 5° -^ {pirf decays [38]. 

We have also presented an improved measurement of the CP-violating charge asymmetry Ak-k 
in the B^ -^ K'^tt" decay. We observe direct CP violation in B^ -^ K~^7r~ with a significance 
of 6.I0". Ignoring color-suppressed tree amplitudes, the charge asymmetries in K~^7:~ and K'^tt^ 
should be equal (see Gronau and Rosner in Ref. [21]), which has not been supported by recent 
BABAR and Belle data [5,7,13]. These results might indicate a large color-suppressed amplitude, 
an enhanced electroweak penguin, or possibly new-physics effects [39]. 

Finally, we have presented an improved measurement oiB{B^ -^ K^t:^). From the rate sum-rule 
prediction [21] 2B{K'^tt°Y' = B{K+tt-) + ^[B{K'^tt+) - 2B{K+tt^)] and the currently published 
results for the other three B -^ Ktt modes, we find the sum-rule prediction to be B{B^ -^ K^ir^y^ = 
(8.4 ± 0.8) X 10^^, which is consistent with our new experimental result. 
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Figure 9: (Top) constraint on the angle Aotttt = a— Oeff, expressed as one minus the confidence level 
(C.L.), as a function of | AoTTTrl- We find an upper bound on |Aa7r7r| of 43° at the 90% C.L. (Bottom) 
constraint on the CKM angle a expressed as 1— C.L. There are eight peaks, two of them nearly 
merged, corresponding to an eight-fold ambiguity in the extraction of a; four solutions are from 
the value and sign of Aaj^T^, which is doubled due to the trigonometric reflections between Oes and 
7r/2 — Oes- We exclude the range [23°, 67°] in a at the 90% C.L. Only the isospin-triangle relations 
and the expressions in Eq. 1 are used in this constraint. The point a = 0, which corresponds to no 
CP violation, and the values of a near or tt can be excluded with additional physics input [13,36]. 
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